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Изложен новый метод решения систем линейных алгебраических уравнений, основанный на применении 
теории меры в вероятностном пространстве; показана возможность его применения для решения уравнения 
линейного межотраслевого баланса. 
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Введение. Задача оптимального планирования является одной из фундаментальных задач со- 
временной экономики. Производственные ресурсы могут сочетаться различными способами и 
обеспечивать при этом одинаковое количество выпускаемой продукции. Например, можно произ- 
водить одно и то же количество продукции при больших финансовых затратах, но малых трудо- 
затратах. Либо, наоборот, делать малые капиталовложения при больших затратах труда. Суть 
оптимального планирования состоит в выборе наиболее рационального способа производства 
продукции, который будет направлен на минимизацию затрат и максимизацию прибыли. 

Одним из способов планирования производства является модель многоотраслевой эконо- 
мики Леонтьева, а точнее, уравнение межотраслевого баланса. С его помощью удобно планиро- 
вать валовой выпуск продукции для какого-либо периода времени по заданному вектору конечно- 
го потребления. Уравнение межотраслевого баланса представляет собой систему линейных ал- 
гебраических уравнений и решается такими методами, как метод Гаусса, Жордана — Гаусса, Кра- 
мера, Якоби, матричный и итерационный методы и др., каждый из которых имеет свои достоинст- 
ва и недостатки. 

Целью данного исследования является разработка нового стохастического метода реше- 
ния уравнения межотраслевого баланса, основанного на использовании теории меры в вероятно- 
стном пространстве и метода простых итераций. 

Суть модели Леонтьева и методы решения уравнения межотраслевого баланса. Модель 
Леонтьева используется для расчёта связи между отраслями через выпуск и потребление продук- 
ции разного вида [1]. На первом этапе строятся уравнения, описывающие балансовые соотноше- 
ния между отраслями (1). Балансовый принцип связи различных отраслей промышленности со- 
стоит в том, что валовой выпуск Ёй отрасли должен быть равным сумме объёмов потребления в 
производственной и непроизводственной сферах. В самой простой форме балансовые соотноше- 
ния имеют вид: 
Ху =Хи+ХЬ +. + Хи +У,; = 2,...П, (1) 
где х; — валовой выпуск Ёй отрасли; х, — объём продукции Ёй отрасли, потребляемый }й отрас- 
лью при производстве объёма продукции х; у, — объём продукции Ёй отрасли, предназначен- 
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ный для реализации (потребления) в непроизводственной сфере, или так называемый про- 
дукт конечного потребления. 
В соотношениях баланса принято рассматривать стоимостный баланс, поскольку продук- 
ция различных отраслей имеет разные измерения. 
На основе анализа экономики США В. В. Леонтьевым был установлен факт, что в течение 


Хх, 

длительного времени соотношения —^ очень слабо меняются и мо! ут рассматриваться как посто- 
Х; 
й 


янные величины. 
В силу указанного факта можно сделать следующее допущение: для производства про- 
дукции /-Й отрасли объёма х; нужно использовать продукцию Ёй отрасли объёма а;х, где 


=" — постоянное число. При таком допущении технология производства принимается ли- 
й 

нейной, а само допущение называется гипотезой линейности. При этом числа а; называются ко- 

эффициентами прямых затрат. Согласно гипотезе линейности имеем: 


х, =а,х,; 1=12,...,П. (2) 


Й 


Тогда уравнения балансовых соотношений можно записать в матричной форме: 


х=Ахч+у, (3) 
где 
Хх, у, а, а,› а,„ 
5-Х — У а» а а 
2 2 21 22 2 
Х = ‚Уу= ‚А= Е (4) 
Х» У» а, а, их Я 


Описание матричного представления (4) и уравнение (3), собственно, и представляют со- 
бой модель Леонтьева. А уравнение (3) носит название уравнения линейного межотраслевого ба- 
ланса. 

Для его решения используются традиционные методы решения систем линейных алгеб- 
раических уравнений (СЛАУ), такие как матричный метод, метод Гаусса, Жордана — Гаусса, Кра- 
мера, простых итераций и т. д. Они подробно изложены во многих источниках, в частности в [2, 
3]. Авторами данной работы предлагается использовать для решения уравнения (3) новый метод, 
имеющий ряд достоинств по сравнению с ранее известными, как применительно к решению урав- 
нения линейного межотраслевого баланса, так и к другим задачам, в которых решаются СЛАУ. 
Применение стохастического метода решения СЛАУ к уравнению межотраслевого ба- 
ланса. Рассмотрим линейное уравнение: 


2(х)=С(х)+у, (5) 
где хеО (С) — искомый входной элемент; 2 (х) =Е — выходной элемент; С— линейный опера- 
тор, в частности для конечномерного пространства, С — матрица пхп элементов; уеЁ — 
вектор правых частей. 
Линейный оператор Сопределён на линейном многообразии (С ) банахова пространства 
Еи отображает его в некое банахово пространство 2. 
Множество 2 (С) — область определения оператора С. Совокупность всех 2 = , для ко- 


торых уравнение (5) разрешимо, является областью значений оператора С. 
В уравнении (5) мы сталкиваемся с отображением, которое в конечномерном пространстве 
сужается до векторной функции случайных аргументов. Необходимо определить вероятностные 
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характеристики векторной случайной величины 2 по заданным вероятностным характеристикам 
векторной случайной величины Х, связанной с величиной 2 функциональной зависимостью 
2 = Ф(Х) : 

Для того чтобы полностью охарактеризовать случайную величину, необходимо задать 
множество её возможных значений и дать способ определения и сравнения между собой вероят- 
ностей этих значений. Такая характеристика случайной величины называется её законом распре- 
деления. 

Наиболее общей формой закона распределения случайной величины является её функция 
распределения [4]. Она имеет вид: 

Р (ха ЕР (О (5х, 

Функция Е (х) является вероятностной мерой. 


Если 2 — функция векторного случайного аргумента, то её функция распределения имеет 
вид: 
Е (2)=Р(2 32) ЕР (Хе) = [д | Рон)боасаю (6) 

(хо)  Ф(Х»)< 2 
где интеграл является п-кратным и распространён на все области, для которых выполнено нера- 
венство, написанное под знаком интеграла. 

В частном случае, когда функция ф монотонна в области возможных значений случайной 
величины Х, существует только одна область, в которой выполнимо неравенство под знаком ин- 
теграла (6). В этом случае интеграл заменой переменных приведётся к виду: 


г (2) 5 | | Г (ч (о и (п) ,-- м (Пр ))- 7 апобй,...а, , (7) 


Ох 


2ч(п)|_|2х 
р2 


р 


Если матрица Якоби не обращается в ноль в области возможных значений, то существует 


где п=Фх(х); х=у (п); /— матрица Якоби вида 1 = 














„0х (02) 
обратное отображение, причём — =| — |. 
22 рх 
Дифференцируя (7), получим формулу, характеризующую закон сохранения меры для 
плотности вероятности векторной случайной величины 2: 


Г, (2) =" (у(2)). (8) 
В исследуемом авторами случае векторная величина 2 является линейной функцией век- 
торной случайной величины Х, поэтому 


2=Ф(х)=Ск+у, х=у(2)=С".2-С“.у. (9) 


Подставляя последнее выражение в (8), получим формулу для плотности вероятности ли- 
нейной функции векторной случайной величины: 


Е, (2) =И(С*.2-с7 у). (10) 


Эта формула показывает, что при линейном преобразовании векторной случайной вели- 
чины кривая плотности вероятности не изменяет своего характера, а только смещается относи- 
тельно начала координат. 

Применим формулу (10) для решения уравнения (5). 

Как было сказано выше, предлагаемый в статье подход использует метод итераций, по- 
этому необходимо обозначить начальную область изоляции корня [% и привести уравнение (5) к 
виду: 
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х=Ахчу, (11) 
где А=С+Е, Е— единичная матрица. 
Как видно, уравнение (11) представляет собой не что иное, как уравнение линейного 
межотраслевого баланса (3). 
Для сходимости итерационного процесса достаточно, чтобы оператор А из уравнения (11) 
был сжимающим [5, 6], т. е. |А|<1. 


Если оператор А не сжимающий, то его необходимо привести к сжимающему виду. 
Каждое последующее приближение х будем находить по формуле 


х, =А- Хх, +У. (12) 
Из (5) следует, что 2(х„.)=С.х„, +У , поэтому 
Хх, =А- хачу =А?.х,,+(А+Е)-у=А-х, + (А +А+Е)-у=... 
Таким образом, можно выразить приближение хна любой итерации через х формулой 
х, = А". Жо + (А+ А"? +...+Е).у. (13) 


Обозначим (4 но Е). у=р. Тогда 


Хх, =А".х, +, х=А"-(х,„-0). (14) 


Подставим последнее уравнение в формулу (10): 


| о" -2}} _ г (= 


п 





сои" 








п р (15) 
здесь внешние прямые скобки в определителе означают, что этот определитель берётся по абсо- 
лютной величине. 

С учётом преобразований 


Е, (х,) = Е, (А. (х, - А” - А"? -...-Е)-у)-|А" 


Хя < п 








р (16) 





где Г, (^” :(х ДА иЕЕЕ ). у) — значение плотности вероятностной меры начального 


Ш 


приближения ж в области изоляции корня [%. 

Область [% задаётся 2п точками, у каждой из которых л координат. Координаты точек об- 
ласти О на последующих итерациях определяются по формуле (12). 

Область Р сжимается к точке концентрации корня и вероятностная мера вырождается в 
этой точке в 5-функцию, если считать оператор и вектор правых частей не случайными. Каждое 
найденное по формуле (15) значение плотности вероятностной меры на Ёй итерации имеет место 
в области 0. 

Из полученной в результате подстановки формулы (16) видно, что, задав начальное при- 
ближение и определив плотность его вероятностной меры, можно найти плотность вероятностной 
меры хна любом шаге. 

Заключение. В работе изложен принципиально новый метод решения систем линейных алгеб- 
раических уравнений и рассмотрено его применение для решения уравнения линейного межот- 
раслевого баланса. 

Предложенный подход выгодно отличается от используемых ранее методов тем, что в 
анализе участвует достаточно широкая окрестность решения. Также неоспоримым преимущест- 
вом метода является отсутствие затруднений с выбором начальной точки. 
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Кроме того, предложенный алгоритм имеет более широкую область применения, чем ра- 
нее известные. В частности, он может быть использован в случае нечёткого задания линейного 
оператора или вектора правых частей. 
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